Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Seven-parameter statistical model for BRDF in the UV band

Open Access Open Access

Abstract

A new semi-empirical seven-parameter BRDF model is developed in the UV band using experimentally measured data. The model is based on the five-parameter model of Wu and the fourteen-parameter model of Renhorn and Boreman. Surface scatter, bulk scatter and retro-reflection scatter are considered. An optimizing modeling method, the artificial immune network genetic algorithm, is used to fit the BRDF measurement data over a wide range of incident angles. The calculation time and accuracy of the five- and seven-parameter models are compared. After fixing the seven parameters, the model can well describe scattering data in the UV band.

©2012 Optical Society of America

1. Introduction

Interest in the study of scattering and reflective properties in the UV band (200 nm to 400 nm) has recently increased. The applications of such properties include UV space object detection and subsurface defect short wavelength monitoring in the semiconductor industry, among others [1,2].

The bidirectional reflectance distribution function (BRDF) is often used to describe the directional dependence of the scattering properties of a surface. The BRDF has been extensively studied and surveyed in various fields [310]. Many BRDF models have been developed, and they can be classified into two categories, namely, purely empirical and analytical models.

Both models have their advantages and disadvantages. The purely empirical BRDF models are simple and useful but have no physical basis, such as the Minnaert BRDF model [3], Walthall model [5], etc. The analytical BRDF models are derived from more complex physical theory by simplifying some assumptions and approximations. They are needed in many applications, especially in simulation software, but they usually require the determination of many parameters. Examples include the Roujean [7], LiSparse-Dense BRDF [8], and Renhorn and Boreman 14-parameter [9] models.

Most BRDF models are developed over the visible and near-infrared spectrum. There are numerous papers that have discussed the measurement or model method in the visible to infrared band [712]. Although there are some related papers reporting UV scatter in rough surfaces [1317], a suitable BRDF model in the UV band is still lacking. In the UV band, the refraction light spectrum is often not easy to be measured, and we should consider more complex model to describe the measurement data. Bulk or volume scattering are usually caused by the in-homogeneities of materials. Relative studies about volume scattering provide a powerful tool to investigate the scattering properties of roughness surfaces [1822].

The present paper aims to provide a model of reflection by a rough surface that can successfully predict the experimental findings in the UV band. In section one, the schematic diagram of the instrument used to performing scatter measurements is presented. In section two, a novel seven-parameter model is developed. There are three terms in this model: surface scatter, bulk scatter, and retro-reflection. An optimizing modeling method is used to model the BRDF measurement data of typical samples in the UV band. In the last section, the calculation results using the new seven-parameter model are compared with the five-parameter model.

2. Measurement of BRDF in the UV band

A BRDF is actually an angle-resolved energy distribution. As shown in Fig. 1 , an element of surface dA is illuminated by a source of incident wave with wave vectorki. Letting kr be the reflected direction wave vector, the symbol z^ denotes the normal of the mean surface of the z axis and n^ denotes the normal direction of the micro-facetdA. Letting α be the angle between n^ and z^, where γ is the angle between ki and n^, angles α and γ satisfy the following relationships

cosα=(cosθi+cosθr)/(2cosγ)
cos2γ=12(cosθicosθr+sinθisinθrcosφr+1)
where θ and φ are the zenith and azimuthal angles, respectively. The subscripts i and r denote the incident and reflected angles, respectively. The definition of BRDF,fr, can be expressed as the differential radiance dLr scattered by a uniformly illuminated, homogeneous material per unit differential incident irradiance dEi

 figure: Fig. 1

Fig. 1 Angle definition of BRDF.

Download Full Size | PDF

fr(θi,φi,θr,φr)=dLr(θi,φi,θr,φr)dEi(θi,φi)

Another equivalent definition of the BRDF is

fr=Ps/(PicosθrΩs)
where Ωs is the solid angle through which the scattered power Ps is collected. It is then normalized with respect to the total incident powerPi. On the other hand, cosθr is the factor that can be thought of as a correction to adjust the illuminated area to its apparent area when viewed from the direction of the scatter.

Figure 2 shows the schematic diagram of the instrument used to perform angle-resolved optical scatter measurements. This BRDF measurement installation system was designed and constructed by the Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences. This measurement system has three parts, i.e., the source, sample holder, and receiver. The design details of each part are described elsewhere [23].

 figure: Fig. 2

Fig. 2 A schematic diagram of the instrument.

Download Full Size | PDF

The UV light source used in these investigations was a 150 mW-power deuterium lamp, which can generate continuous emission spectra. After emerging from the exit slit, the beam is quasi-collimated by passing through a convex lens. After passing through two mirrors, the beam is reflected to the slit. By passing through the slit, the beam enters the sample and is then be collected by the detector. The system can measure three-dimensional BRDF data by rotating the sample and light source using electromotors A, B, and C in different three-dimensional directions. The detector used was a spectrometer made by Ocean Optic Co. in USA. The sample can be rotated about a vertical axis to allow different angles of incidence, and the detector can be rotated independently about the same axis to allow the measurement of angular scatter distribution. The step of the zenith angle is 1°, and the step of the azimuth angle is 5°. Near the specular direction, the sample angle is 1° to guarantee sufficient precision. The accuracy of the zenith and azimuth angles is 0.1°, and the relative error of the measured BRDF of this system is less than 5%. The incident UV wavelength varies from 250.45288 nm to 368.89652 nm. The total 1024 output voltages Vs for the sample are measured within this wavelength range. The voltages of reflectance standard plate Vref are also measured within the corresponding wavelengths. Then, the ratio is as follows

frs(θi,φi;θr,φr,λ)frr(θi,φi;θr,φr,λ)=VsVref
where frs(θi,φi;θr,φr,λ) and frr(θi,φi;θr,φr,λ) are the BRDFs of the sample and reflectance standard plate at a given incident wavelengths λ, respectively.

Notably, Vref of the reference standard plate must be determined under the same measurement conditions used to detect the samples, and all measurements must be performed correspondingly. In our experiments, the reference standard plate is the pressed polytetrafluoroethylene plate, which can be considered as a perfect Lambertian plate.

3. Development of a seven-parameter BRDF model

Torrance and Sparrow [10] have deduced a useful light-scattering model, which assumes that the surface consists of small, randomly disposed mirror-like facets. An implied criterion of this model is that the root-mean-square (RMS) surface roughness is greater than the wavelength of the incident radiant energy. There are three parameters in the Torrance-Sparrow model, based on which Wu [11] has proposed a five-parameter BRDF model, whose expression can be expressed as

BRDFsurface=fr(θi,θr,φ,λ)=kbkr2cosα1+(kr21)cosαexp[b(1cosγ)a]V(θi,θr,φ,λ)cosθicosθr+kdcosθi
where φ=φrφi=φr. For isotropic materials at a given wavelength, the BRDF can be safely assumed to be invariant to azimuthal rotations of the incident andφi=0. There are two components in Eq. (6); one is the specular reflection from mirror-like surface facets, and the other is a diffuse component. In the modified Torrance-Sparrow model of Wu [11], there are five parametersa, b, kb, kd, and kr to be determined. kbis the mirror-direction component, kd is related with the diffuse reflection component, and kr is related to the distribution of the subsurface dA (determined by the slope distribution of the surface). kr2cosα/[1+(kr21)cosα]is the distribution function of the sub-surface, exp[b(1cosγ)a] is the approximation description of the Fresnel reflectance function, and V(θi,θr,φ,λ) is the masking and shadowing effect. In this model, the exponential function in the Torrance-Sparrow model is substituted by an elliptical function to describe the distribution of the normal of the facets. The exponential function with two parameters is used to substitute the Fresnel reflectance function to avoid the calculation of many trigonometric functions. The model can be used to describe isotropic surfaces with non-polarization incident light. The selection criteria of these five parameters are to minimize the RMS errors between the simulation and experimental data.

UV wavelengths are shorter than visible and infrared bands. Thus, a new model that suited the experimental measurement data in the UV band is developed. Our model partly adopted the five-parameter model of Wu and partly adopted the Renhorn-Boreman fourteen-parameter analytical model. The corresponding bulk scatter and retro-reflection terms are introduced to improve the accuracy of the simulation results of the five-parameter model.

According to the Renhorn-Boreman model [9], a Gaussian surface auto-covariance function g(x,y) is given by

g(x,y)=σg2exp[(x2+y2)ρ2]
where σg is the RMS surface roughness and ρ is the inverse of the surface correlation length. Its Fourier transform corresponding to the surface power spectrum isG(ξ,η). The parameters ξ and η are direction cosines defined asξ0=sin(θi), ξ=sin(θr)cos(φ) andη=sin(θr)sin(φ). The corresponding Fourier transform of the surface auto-covariance function is

G(ξ,η)=σg22ρexp[14ρ2(ξ2+η2)]

An exponential auto-covariance function is given by

g(x,y)=σg2exp[(|x|+|y|)ρ]

The Fourier transformation results of Eq. (9) can be obtained as follows using Mathematic under the condition 0ξ1,0η1 and a sufficiently largeρ

G(ξ,η)=2πσg2ξ2+η2+ρ2

In the Renhorn-Boreman model, the exponential statistical distribution is combined with the Gaussian distribution. Taking the angle of incidence into account, the Renhorn-Boreman model yields the two-dimensional Lorentzian BRDF by

BRDF=σN(ξξ0)2+η2+ρ2
where σ is the integrated reflectance. To ensure that the analytical two-dimensional integral of the BRDF is unitary, the parameter N gives

N=1πlog[ξ02+ρ2+ρ4+2ξ02ρ2+2ρ2+(1ξ02)2+12ρ2]

Considering the geometrical foreshortening effect, the expression ρ(1ξ02+1ξ2η2)/2 is used instead of the surface correlation parameterρ. Equation (11) of the BRDF model of Renhorn-Boreman is finally given as

BRDF=σN(ξξ0)2+η2+14(1ξ02+1ξ2η2)2ρ2

The Renhorn-Boreman BRDF model, which considers the surface scattering and includes the bulk and retro-reflection scattering sections, is partly adopted. The total BRDF is given by three separate parts as follows

BRDFtotal=BRDFsurface+BRDFbulk+BRDFretro
The surface scattering BRDF part is described by Eq. (6), whereas the bulk scattering and retro-reflection scattering BRDF are described as follows
BRDFbulk=N(ξξ0)2+η2+14(1ξ02+1ξ2η2)2ρ12
BRDFretro=N(ξξ0)2+η2+14(1ξ02+1ξ2η2)2ρ22
where ρ1 and ρ2 are the parameters related to the bulk and retro-reflection scatterings, respectively. Given that the proportions of the bulk and retro-reflection scattering sections are very small, only the masking and shadowing effects are considered in the surface scattering section. For the bulk and retro-reflection scattering sections, only the two-dimensional Lorentzian BRDF is considered.

4. Comparison of the seven- and five-parameter BRDF models

At this point, a new BRDF model has been developed. There are seven parameters to be determined in this model, namely, {kb,kd,kr,a,b,ρ1,ρ2}. The artificial immune network algorithm is used to fit the parameters in the model. The optimization method is used to determine the optimum values of the fitting parameters of the BRDF to minimize the squared error between the measured data and the model.

  • (1) Measured BRDF data of sample #1: Fine grinding rough aluminum plate.

Figures 3 and 4 compare the measured data (shown as discrete points) with the fitted BRDF derived from the seven-parameter (solid curves in blue) and five-parameter (short dash curves in red) models. The sample used in this simulation is a fine grinding rough aluminum plate. The zenith angles of incidence are 20°, 30°, 45°, and 60°. The scattering angle of scatter is over a range of angles from –70° to 70°. In the given wavelengths, six groups of measured data at various angles of incidence are usually used to fit the parameters to be determined (only four groups of measured data are shown in the Figs. 3-6 ). For convenience, the azimuthal angles of incidence and scatter are both set to 0°. The wavelengths of incidence are λ=266nm and λ=369nm in Figs. 3 and 4, respectively.

 figure: Fig. 3

Fig. 3 Comparison of the seven- and five-parameter models of sample #1 for various angles of incidence when wavelength of incidence isλ=266nm.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 Comparison of the seven- and five-parameter models of sample #1 for various angles of incidence when wavelength of incidence isλ=369nm.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 Comparison of the seven- and five-parameter models of sample 2# for various angles of incidence when wavelength of incidence isλ=266nm.

Download Full Size | PDF

 figure: Fig. 6

Fig. 6 Comparison of the seven- and five-parameter models of sample 2# for various angles of incidence when wavelength of incidence isλ=369nm.

Download Full Size | PDF

Based on Figs. 3-4, sample #1 can be concluded as sufficiently rough, and the diffuse component dominates this rough surface. Compared with the same sample at different incident wavelengths, Figs. 3 and 4 show that if the wavelength is shorter, the sample is rougher. The mirror reflection at 369 nm is distinct from that at 266 nm. At both wavelengths, the seven-parameter model can meet the measured data better than the five-parameter one.

Table 1 lists the fitting parameters, RMS error, and calculation times of the two different models. Compared with the five-parameter model, our newly developed seven-parameter model has a smaller RMS error. Although the calculation time is longer than that of the five-parameter model, it is still tolerable.

Tables Icon

Table 1. Comparison of the Five- and Seven-Parameter Models of Sample #1

All calculation time simulations shown in Table 1 is developed using a Fortran code running on a standard computer (Intel Core 2 Duo processor, 1.86 GHz, 500GB EMS memory). The parameters are determined by an optimizing artificial immune network genetic algorithm.

According to the Kirchhoff approximation method for electromagnetic wave scattering, the bulk scattering increases with decreased incident wavelength. This model considers bulk and retro-reflection scatterings. Thus, it can effectively reduce the RMS error between the experimental data and model simulation results.

(2) Measured BRDF data of sample #2: white painted surface

Another example is shown in Figs. 5 and 6. Table 2 lists the model parameters, RMS error, and computing time of sample #2. The sample used in this case is a white painted surface.

Tables Icon

Table 2. Comparison of the Five- And Seven-Parameter Models of Sample #2

The mirror-reflective section of sample #2 is more obvious than that of sample #1. Thus, sample #2 is relatively smoother than sample #1. The seven-parameter model still has less RMS error for different incidence wavelengths in this case.

5. Conclusion

A new seven-parameter BRDF model of rough surfaces in the UV band is established. An angle-resolved BRDF measurement system operates on some typical samples at different incident angles. The measured data are compared with the results of the seven-parameter model. Both for different wavelengths of the same sample or for different samples, the new seven-parameter model meets the experimental data better than the five-parameter one. It can significantly improve the calculation accuracy, although it costs a little more calculation time.

Both seven- and five-parameter model has its advantages and disadvantages. Optimization method will become more difficult when parameter in a model is so many. And the calculation time will be longer. If the precision requirement is not high, for most engineering-oriented application fields, the five-parameter model is a good enough model. But for some special use such as in the UV band, it is difficult to detect the scattering signal. If you want to improve the calculation precision, more precision model should be used.

This paper offers a new developed BRDF model of rough surfaces in the UV band. The scattering properties of materials in the UV band are elucidated. And this kind of researches is helpful to understand the scattering properties of materials in the UV band.

Acknowledgments

This work was funded by the National Natural Science Foundation of China (No. 60971065, 61172031) and the Fundamental Research Funds for the Central Universities.

References and links

1. C. Lavigne, G. Durand, and A. Roblin, “Ultraviolet light propagation under low visibility atmospheric conditions and its application to aircraft landing aid,” Appl. Opt. 45(36), 9140–9150 (2006). [CrossRef]   [PubMed]  

2. M. A. Velazco-Roa and S. N. Thennadil, “Estimation of complex refractive index of polydisperse particulate systems from multiple-scattered ultraviolet-visible-near-infrared measurements,” Appl. Opt. 46(18), 3730–3735 (2007). [CrossRef]   [PubMed]  

3. M. Minnaert, “The reciprocity principle of linear photometry,” Astrophys. J. 93, 403–410 (1941). [CrossRef]  

4. J. Stover, Optical Scattering, Measurement and Analysis (SPIE Press, 1995).

5. C. L. Walthall, J. M. Norman, J. M. Welles, G. Campbell, and B. L. Blad, “Simple equation to approximate the bidirectional reflectance from vegetative canopies and bare soil surfaces,” Appl. Opt. 24(3), 383–387 (1985). [CrossRef]   [PubMed]  

6. P. Beckman and A. Spizzichino, The Scattering of Electromagnetic Waves from Rough Surfaces (Pergamon, 1963).

7. J. L. Roujean, M. Leroy, and P. Y. Deschamps, “A bidirectional reflectance model of the earth’s surface for the correction of remote sensing data,” J. Geophys. Res. 97(20), 455–468 (1992).

8. X. Li and A. H. Strahler, “Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: effect of crown shape and mutual shadowing,” IEEE Trans. Geosci. Rem. Sens. 30(2), 276–292 (1992). [CrossRef]  

9. I. G. E. Renhorn and G. D. Boreman, “Analytical fitting model for rough-surface BRDF,” Opt. Express 16(17), 12892–12898 (2008). [CrossRef]   [PubMed]  

10. K. E. Torrance and E. M. Sparrow, “Theory for off-specular reflection from roughened surfaces,” J. Opt. Soc. Am. 57(9), 1105–1114 (1967). [CrossRef]  

11. Z. S. Wu, D. H. Xie, P. H. Xie, and Q. N. Wei, “Modeling reflectance function from rough surface and algorithms,” Acta Opt. Sin. 22, 897–901 (2002).

12. K. J. Dana, B. van Ginneken, S. K. Nayar, and J. J. Koenderink, “Reflectance and texture of real-world surfaces,” ACM Trans. Graph. 18(1), 1–34 (1999). [CrossRef]  

13. Y. Barnes and J. J. Hsia, “UV bidirectional reflectance distribution function measurements for diffusers,” Proc. SPIE 1764, 285–288 (1993). [CrossRef]  

14. M. P. Newell, L. A. Whitlock, and R. A. M. Keski-Kuha, “Extreme ultraviolet scatter from particulate contaminated mirrors,” Proc. SPIE 2541, 174–185 (1995). [CrossRef]  

15. M. P. Newell and R. A. M. Keski-Kuha, “Bidirectional reflectance distribution function of diffuse extreme ultraviolet scatterers and extreme ultraviolet baffle materials,” Appl. Opt. 36(22), 5471–5475 (1997). [CrossRef]   [PubMed]  

16. T. Zurbuchen, P. Bochsler, and F. Scholze, “Reflection of ultraviolet light at 121.6 nm from rough surfaces,” Opt. Eng. 34(5), 1303–1315 (1995). [CrossRef]  

17. M. P. Newell and R. A. M. Keski-Kuha, “Extreme ultraviolet BRDF measurements: instrumentation and results,” Proc. SPIE 2864, 453–464 (1996). [CrossRef]  

18. C. Amra, “Light scattering from multilayer optics. I. Tools of investigation,” J. Opt. Soc. Am. A 11(1), 197–210 (1994). [CrossRef]  

19. C. Amra, “Light scattering from multilayer optics. II. Application to experiment,” J. Opt. Soc. Am. A 11(1), 211–226 (1994). [CrossRef]  

20. C. Amra, “From light scattering to the microstructure of thin-film multilayers,” Appl. Opt. 32(28), 5481–5491 (1993). [CrossRef]   [PubMed]  

21. J. M. Elson, J. P. Rahn, and J. M. Bennett, “Light scattering from multilayer optics: comparison of theory and experiment,” Appl. Opt. 19(5), 669–679 (1980). [CrossRef]   [PubMed]  

22. C. Amra, C. Grèzes-Besset, and L. Bruel, “Comparison of surface and bulk scattering in optical multilayers,” Appl. Opt. 32(28), 5492–5503 (1993). [CrossRef]   [PubMed]  

23. H. L. Zhang, Z. S. Wu, Y. H. Cao, and G. Zhang, “Measurement and statistical modeling of BRDF of various samples,” Opt. Appl. 40, 197–208 (2010).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 Angle definition of BRDF.
Fig. 2
Fig. 2 A schematic diagram of the instrument.
Fig. 3
Fig. 3 Comparison of the seven- and five-parameter models of sample #1 for various angles of incidence when wavelength of incidence is λ=266nm .
Fig. 4
Fig. 4 Comparison of the seven- and five-parameter models of sample #1 for various angles of incidence when wavelength of incidence is λ=369nm .
Fig. 5
Fig. 5 Comparison of the seven- and five-parameter models of sample 2# for various angles of incidence when wavelength of incidence is λ=266nm .
Fig. 6
Fig. 6 Comparison of the seven- and five-parameter models of sample 2# for various angles of incidence when wavelength of incidence is λ=369nm .

Tables (2)

Tables Icon

Table 1 Comparison of the Five- and Seven-Parameter Models of Sample #1

Tables Icon

Table 2 Comparison of the Five- And Seven-Parameter Models of Sample #2

Equations (16)

Equations on this page are rendered with MathJax. Learn more.

cosα= ( cos θ i +cos θ r ) / (2cosγ)
cos 2 γ= 1 2 ( cos θ i cos θ r +sin θ i sin θ r cos φ r +1 )
f r ( θ i , φ i , θ r , φ r )= d L r ( θ i , φ i , θ r , φ r ) d E i ( θ i , φ i )
f r = P s / ( P i cos θ r Ω s )
f rs ( θ i , φ i ; θ r , φ r ,λ) f rr ( θ i , φ i ; θ r , φ r ,λ) = V s V ref
BRD F surface = f r ( θ i , θ r ,φ,λ) = k b k r 2 cosα 1+( k r 2 1 )cosα exp[ b ( 1cosγ ) a ] V( θ i , θ r ,φ,λ ) cos θ i cos θ r + k d cos θ i
g(x,y)= σ g 2 exp[ ( x 2 + y 2 ) ρ 2 ]
G(ξ,η)= σ g 2 2 ρ exp[ 1 4 ρ 2 ( ξ 2 + η 2 ) ]
g(x,y)= σ g 2 exp[ (| x |+| y |)ρ ]
G(ξ,η)= 2 π σ g 2 ξ 2 + η 2 + ρ 2
BRDF= σN ( ξ ξ 0 ) 2 + η 2 + ρ 2
N= 1 πlog[ ξ 0 2 + ρ 2 + ρ 4 +2 ξ 0 2 ρ 2 +2 ρ 2 + (1 ξ 0 2 ) 2 +1 2 ρ 2 ]
BRDF= σN (ξ ξ 0 ) 2 + η 2 + 1 4 ( 1 ξ 0 2 + 1 ξ 2 η 2 ) 2 ρ 2
BRD F total =BRD F surface +BRD F bulk +BRD F retro
BRD F bulk = N (ξ ξ 0 ) 2 + η 2 + 1 4 ( 1 ξ 0 2 + 1 ξ 2 η 2 ) 2 ρ 1 2
BRD F retro = N (ξ ξ 0 ) 2 + η 2 + 1 4 ( 1 ξ 0 2 + 1 ξ 2 η 2 ) 2 ρ 2 2
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.